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Abstract 

We calculate the cross section of J/ip plus jet inclusive production in 77 colli- 
sions at next-to-leading order within the factorization formalism of nonrelativistic 
quantum chromodynamics (NRQCD) focusing on direct photoproduction. Apart 
from direct J/ip production, we also include the feed-down from directly-produced 
Xcj and ip' mesons. We discuss the analytical calculation, in particular the treat- 
ment of the various types of singularities and the NRQCD operator renormalization, 
in some detail. We present theoretical predictions for the future e + e~ linear collider 
TESLA, taking into account both brems- and beamstrahlung. 

PACS numbers: 12.38.Bx, 12.39. St, 13.66.Bc, 14.40. Gx 



1 Introduction 



Since the discovery of the J/if> meson in 1974, charmonium has provided a useful labora- 
tory for quantitative tests of quantum chromodynamics (QCD) and, in particular, of the 
interplay of perturbative and nonperturbative phenomena. The factorization formalism 
of nonrelativistic QCD (NRQCD) jl!2j provides a rigorous theoretical framework for the 
description of heavy- quarkonium production and decay. This formalism implies a separa- 
tion of short-distance coefficients, which can be calculated perturbatively as expansions in 
the strong-coupling constant a s , from long-distance matrix elements (MEs), which must 
be extracted from experiment. The relative importance of the latter can be estimated by 
means of velocity scaling rules; i.e., the MEs are predicted to scale with a definite power of 
the heavy-quark (Q) velocity v in the limit d<1. In this way, the theoretical predictions 
are organized as double expansions in a s and v. A crucial feature of this formalism is that 
it takes into account the complete structure of the QQ Fock space, which is spanned by 
the states n = 2S + 1 l l y' with definite spin S, orbital angular momentum L, total angular 
momentum J, and colour multiplicity a = 1, 8. In particular, this formalism predicts 
the existence of colour-octet (CO) processes in nature. This means that QQ pairs are 
produced at short distances in CO states and subsequently evolve into physical, colour- 
singlet (CS) quarkonia by the nonperturbative emission of soft gluons. In the limit v — > 0, 
the traditional CS model (CSM) [3| is recovered. The greatest triumph of this formalism 
was that it was able to correctly describe |3] the cross section of inclusive charmonium 
hadroproduction measured in pp collisions at the Fermilab Tevatron [3] , which had turned 
out to be more than one order of magnitude in excess of the theoretical prediction based 
on the CSM. 

Apart from this phenomenological drawback, the CSM also suffers from severe concep- 
tual problems indicating that it is incomplete. These include the presence of logarithmic 
infrared singularities in the 0(a s ) corrections to P-wave decays to light hadrons and in 
the relativistic corrections to S"-wave annihilation [H] , and the lack of a general argument 
for its validity in higher orders of perturbation theory. While the /^-factorization [Zj 
and hard-comover-scattering jS] approaches manage to bring the CSM prediction much 
closer to the Tevatron data, they do not cure the conceptual defects of the CSM. The 
colour evaporation model which is intuitive and useful for qualitative studies, also 
significantly improves the description of the Tevatron data as compared to the CSM [10 . 
However, it does not account for the process-specific weights of the CS and CO contribu- 
tions, but rather assumes a fixed ratio of 1 : 7. In this sense, a coequal alternative to the 
NRQCD factorization formalism is presently not available. 

In order to convincingly establish the phenomenological significance of the CO pro- 
cesses, it is indispensable to identify them in other kinds of high-energy experiments as 
well. Studies of charmonium production in ep photoproduction [TTj, ep |12|13j and uN 
|14|15j deep- inelastic scattering (DIS), e + e~ annihilation in the continuum [TS], Z-boson 
decays [T7j, 77 collisions |18|19|2(J|21| . and 6-hadron decays |22] m ay be found in the 
literature; for reviews, see Ref. Furthermore, the polarization of ip' mesons pro- 

duced directly |2H an d of J/ip mesons produced promptly [25 26], i.e., either directly or 
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via the feed-down from heavier charmonia, which also provides a sensitive probe of CO 
processes, was investigated. Until recently, none of these studies was able to prove or 
disprove the NRQCD factorization hypothesis. However, HI data of ep — > e + J/tp + X in 
DIS at the DESY Hadron Electron Ring Accelerator (HERA) [21! and DELPHI data of 
77 -> J /if; + X at the CERN Large Electron Positron Collider (LEP2) [28J provide first 
independent evidence for it |13|21| . 

The verification of the NRQCD factorization hypothesis is presently hampered both 
from the theoretical and experimental sides. On the one hand, the theoretical predic- 
tions to be compared with existing experimental data are, apart from very few exceptions 
[29 3 0p31| . of lowest order (LO) and thus suffer from considerable uncertainties, mostly 
from the dependences on the renormalization and factorization scales and from the lack 
of information on the nonperturbative MEs. On the other hand, the experimental errors 
are still rather sizeable. The latter will be dramatically reduced with the upgrades of 
HERA (HERA II) and the Tevatron (Run II) and with the advent of CERN LHC and 
hopefully a future e + e~ linear collider (LC) such as the TeV-Energy Superconducting 
Linear Accelerator (TESLA), which is presently being designed and planned at DESY. 
On the theoretical side, it is necessary to calculate the next-to-leading-order (NLO) cor- 
rections to the hard-scattering cross sections and to include the effective operators which 
are suppressed by higher powers in v. 

In this paper, we take a first step in this comprehensive research programme, by 
studying the inclusive production of J /if} mesons in high-energy 77 collisions. As men- 
tioned above, this process was studied at LEP2 j2H], where the photons originated from 
hard initial-state bremsstrahlung. At high-energy e + e~ LCs, an additional source of hard 
photons is provided by beamstrahlung, the synchrotron radiation emitted by one of the 
colliding bunches in the field of the opposite bunch. The highest possible photon ener- 
gies with large enough luminosity may be achieved by converting the e + e~ LC into a 
77 collider via back-scattering of high-energetic laser light off the electron and positron 
beams. 

In order for a J/if> meson to acquire finite transverse momentum (pr), it must be pro- 
duced together with another particle or a hadron jet (J). From coupling and phase-space 
considerations it is evident that J/tf> plus jet associated production yields the dominant 
contribution. In the following, we thus consider the process 77 — > J /if) + j + X , where 
X denotes the hadronic remnant possibly including a second jet. Here, we take j and X 
to be free of charm assuming that charmed hadrons or charmonia besides the J / if> meson 
would be detectable. The process 77 — > J/if> + 7 + X , where 7 represents a prompt 
photon, will be considered in a forthcoming publication [32J. 

The incoming photons can interact either directly with the quarks participating in 
the hard-scattering process (direct photoproduction) or via their quark and gluon content 
(resolved photoproduction). Thus, the process 77 — > J/if> + j + X receives contributions 
from the direct, single-resolved, and double-resolved channels. All three contributions are 
formally of the same order in the perturbative expansion. This may be understood by 
observing that the parton density functions (PDFs) of the photon have a leading behaviour 
proportional to aln(M 2 /An CD ) oc a/a s , where a is the fine-structure constant, M is the 
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factorization scale, and Aqcd is the asymptotic scale parameter of QCD. In the following, 
we focus our attention on the direct channel. The other channels are left for future work. 

The J/if> mesons can be produced directly; or via radiative or hadronic decays of 
heavier charmonia, such as Xcj and if>' mesons; or via weak decays of b hadrons. The 
respective decay branching fractions are B(xco — * J/ip + 7) = (1.02 ± 0.17)%, B(xd 
J ftp + 7) = (31.6 ± 3.2)%, B( X c2 -> J/ip + 7) = (18.7 ± 2.0)%, B{if)' -> J/ip + X) = 
(55.7 ± 2.6)%, and B(B -> J/ip + X) = (1.15 ± 0.06)% [33]. The b hadrons can be 
detected by looking for displaced decay vertices with dedicated vertex detectors, and the 
J /if) mesons originating from their decays can thus be treated separately. Therefore and 
because of the smallness of B(B — > J/if> + X), J /if) production through 6-hadron decay is 
not considered here. The cross sections of the four residual indirect production channels 
may be approximated by multiplying the direct-production cross sections of the respective 
intermediate charmonia with their decay branching fractions to J /if) mesons. 

To summarize, the goal of the present analysis is to calculate the inclusive cross section 
of 77 — > J/if) + X in direct photoproduction for finite values of at NLO within NRQCD 
allowing for the J /if) meson to be promptly produced. The LO result, also including the 
single- and double-resolved contributions, may be found in Refs. |2U|21j and the references 
cited therein. The leading relativistic correction, which originates from the cc Fock state 
n = 

lp(8) 

and is of 0(v A ) relative to the LO result, was evaluated in Ref. P^j. In that 
paper, also J /if) plus dijet associated production in two-photon collisions was studied for 
direct photoproduction and direct J /if) production. The 2^3 partonic cross sections 
obtained therein constitute a starting point for the evaluation of the real radiative correc- 
tions in our present study. They have to be complemented with the corresponding results 
for direct Xcj production. Furthermore, the kinematic cuts imposed in Ref. |19j to sepa- 
rate the final-state objects have to be removed. This yields soft and collinear singularities, 
which are collectively denoted as infrared (IR) singularities. On top of this, the virtual 
radiative corrections, which arise from 2^2 Feynman diagrams 1 involving closed loops, 
have to be added. They involve IR, ultraviolet (UV), and Coulomb singularities. The 
cancellation of all these singularities is nontrivial and requires the UV renormalization of 
masses, couplings, and wave-functions; the redefinition of NRQCD MEs so as to absorb 
IR and Coulomb singularities; the factorization of initial-state collinear singularities in 
photon PDFs; and the operation of the Kinoshita-Lee-Nauenberg theorem [33] on the 
cancellation of final-state IR singularities. Apart from being of general phenomenological 
relevance, our analysis should thus also be of conceptual interest for the theoretical heavy- 
quarkonium community. After all, this is the first time that the full NLO corrections are 
evaluated for an inclusive 2 — > 2 process within the NRQCD framework. 

This paper is organized as follows. In Section [21 we describe our analytical calculation 
in some detail. Specifically, we discuss the structure of the various types of singularities 
and the mechanisms by which they are removed. Lengthy expressions are relegated to 
the Appendix. In Section El we present our numerical results appropriate for the e + e~ 
mode of TESLA, and discuss their phenomenological implications. Our conclusions are 

1 Here, the cc bound state is considered as one particle. 
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k 


J/^,fp' XcJ 


3 
5 
7 


3p(l) 3o(8) 

lc(8) 3c(8) 3 p(8) 



Table 1: Values of k in the velocity-scaling rule \0 H [n\j oc t> fc for the leading cc Fock 
states n pertinent to H = J/ift, Xcj, ■ 



summarized in Section 0] 



2 Analytic results 



We start this section with a few general remarks. In our analytic calculation, we take 
the colour gauge group to be SU(iV c ) with a generic value of N c , which is put equal to 
3 in our numerical analysis. Colour factors appearing in our formulas include Tp = 1/2, 
C F = (Nt - 1) /(2N C ), C A = N c , and B F = (N? - 4) /(4N C ). We work in the fixed- 
flavour- number scheme, with nj = 3 active quark flavours q — u,d, s, which we treat as 
massless. The charm quark c and antiquark c, with mass m, only appear in the final 
state. We denote the fractional electric charge of quark q by e q . 

The cc Fock states contributing at LO in v are specified for H 
Their MEs satisfy the multiplicity relations 



J/ip,Xcj,^' in Tabled 



3p 



(8) 
J 



QXcj 
QXcJ 



3p(l) 
'3 S (8) 



(2J+1) 

(2J + 1) 
(2J + 1) 



QXcO 
QXcO 



3p(8) 
•M) 



3p(l) 

3 S (sy 



which follow to LO in v from heavy-quark spin symmetry. 

We employ dimensional regularization with d = 4 — 2e space-time dimensions to handle 
the UV and IR singularities, and we introduce a 't Hooft mass \x and a factorization mass 
M as unphysical scales. We formally distinguish between UV and IR poles, which we 
denote as 1/euv an d 1/em, respectively. We apply the projection method of Refs. |31|35j . 
which is equivalent to the ci-dimensional matching procedure of Ref. [SH], in order to ex- 
tract the short-distance coefficients that multiply the MEs. However, in order to conform 
with common standards, we adopt the normalizations of the MEs from Ref. [2] rather 
than from Refs. |3U|31j : i.e., the MEs include spin and colour average factors. 

There is only one partonic subprocess at LO, namely 



7(fci) + 7(^2) -> cc 



(p)+g(h) 



(2) 



which is mediated by the diagrams depicted in Fig. The analogous process with n = 
yields a relativistic correction of 0(v 4 ) and was studied in Ref. ^H]- With the four- 
momentum assignments indicated within the parentheses in Eq. (J2J), we have k\ + ki = 
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p + &3, k\ = k 2 = k\ = 0, and p 2 = 4m 2 . We define the Mandelstam variables as 
s = (ki + k 2 ) 2 , t = (ki — p) 2 , and u = (k 2 —p) 2 , so that s + t + u = 4m 2 . For convenience, 
we also introduce 

Si = s — 4m 2 , ti = t — 4m 2 , u\ = u — Am 2 . (3) 

In the NLO analysis, we need to evaluate the cross section of process (J2J) in d di- 
mensions and retain terms of 0(e) because UV counterterms appear in multiplicative 
renormalization. We have 



da = -dPS 2 (k 1 + k 2 ;p, A; 3 )|T | 2 , (4) 
2s 

where the first factor on the right-hand side stems from the flux, % is the LO transition- 
matrix (T) element of process (0), and it is averaged (summed) over the spin and colour 
states of the incoming (outgoing) particles. Here and in the following, we denote the 
Lorentz-invariant iV-particle phase-space element in d dimensions as 

dPS N (P; Pl , ...,p N )= ^ N -^- d \2n) d 6^ [P - X>) II J^h 5 0* " ^) 6 (f?) > 



i=l 



(5) 

where Pi and are the four-momenta and masses of the final-state particles and P is 
the total four-momentum of the incoming particles. In the case of process (J2J), we have 

1 {4ttu 2 sY 

dPS 2 (h + k 2 ; p, k 3 ) = ^-^tz r — — 6{s + t + u - Am 2 )dt x du u (6) 

87TSl (1 — 6) \ tu J 

where T(x) is Euler's T function, so that 

dao 1 / 47r/i 2 s s 

dtidui 167rs 2 r(l — e) I tu 



5(s + t + u-Am 2 )\%\ 2 . (7) 



Note that the spin average factor of an incoming photon or gluon in d dimensions is 
l/(d - 2). The final result is listed in Eq. (1) of Ref. 02]. 

We now move on to the process e + e~ — ► e + e~H + X, where H denotes a generic char- 
monium state. For later use, we now also include the resolved contributions. We select 
the e + e~ centre-of- momentum (CM) frame and denote the nominal e + e~ energy by VS, 
and the transverse momentum and rapidity of the H meson by pt and y, respectively. 
Invoking the Weizsacker- Williams approximation (WWA) |37j and the factorization the- 
orems of the QCD parton model [HE] and NRQCD the differential cross section of 
e + e~ — > e + e~H + X can be written as 



d 2 a(e + e~ -> e + e~H + X) 
dp\dy 



J dx + f y/e (x + ) J dx_f y/e (x-) J dx a f ah {x ai M) 



a,b,d 

x f iXbhh(Xi , m) u° B w» Mab + d) 



6 



where f 1 / e (x±) is the equivalent number of transverse photons radiated by the initial- 
state positrons and electrons, f a / 7 (x a , M) are the PDFs of the photon, {0 H [n}) are the 
MEs of the H meson, da(ab — > cc[n] + d) are the differential partonic cross sections, 
the integrals are over the longitudinal-momentum fractions of the emitted particles w.r.t. 
the emitting ones, and it is summed over a,b = J,g,q, q and d = g,q,q. With the 
definition / 7 / 7 (x 7 ,M) = 5(1 — x 7 ), Eq. (jHJ) accommodates the direct, single-resolved, 
and double-resolved channels. The Mandelstam variables are given by s = x a XbX + X-S, 



t\ = — x+Xa^fSrriT exp(— y), and u\ = — X-XbVSrnT exp(y), where my = yp^ + Am 2 is 

the transverse mass of the H meson. For a given value of x/S*, the accessible phase space 
is defined by 

S-Am 2 



< p T < 



2VS ' 



, S + 4m 2 

\y\ < arcosh — -= . (9) 

2v Sm T 



2.1 Virtual corrections 



Representative examples of diagrams that generate the virtual corrections to the cross 
section of process (J2J) are shown in Fig. |21 They fall into two classes. The diagrams of the 
first class are obtained by attaching one virtual gluon line in all possible ways to the tree- 
level seed diagrams of Fig. ^ Their contribution is proportional to e 2 c . They include self- 
energy, triangle, box, and pentagon diagrams. Loop insertions in external gluon or charm- 
quark lines are accommodated in the respective wave-function renormalization constants 
and are not displayed in Fig. |21 The self-energy and triangle diagrams are in general 
UV divergent; the triangle, box, and pentagon diagrams are in general IR divergent. 
The pentagon diagrams without three-gluon vertex also contain Coulomb singularities, 
which are cancelled after taking into account the corresponding corrections to the operator 



O 



3 S\ 



(8) 



by a smal 



In the practical calculation, the Coulomb singularities are first regularized 
gluon mass. This regularization prescription is then transformed into one 
implemented with a small relative velocity v between the c and c quarks [29 . 

The diagrams of the second class contain a light-quark triangle or box. The triangle 
diagrams produce contributions proportional to e 2 , or e q e c , which vanish individually by 
Furry's theorem [39J. The box diagrams are proportional to e 2 . They contain UV and IR 
singularities, but their sum is finite. 

Calling the T-matrix element comprising all these loop diagrams %, the virtual cor- 
rections to the cross section of process (J2J) may be evaluated as 

da v 1 / 4-7T /i 2 s' 



dt\dui 167rs 2 r(l — e) y tu 



5 (s + t + u-Am 2 )2Re (%*%). (10) 



We apply two independent approaches to calculate the one-loop diagrams. The first 
one uses FeynArts [IHj to generate the diagrams and self- written Mathematica codes to 
apply the projectors and provide expressions, which are afterwards treated with a FORM 
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program to perform the tensor reduction and the extraction of the UV and IR singularities. 
The result is then transformed into a Fortran code to be used for the numerical evaluation. 
The second approach utilizes QGRAF jH] for the generation of the diagrams, FeynCalc jl2j 
for the tensor reduction, and LoopTools for the numerical evaluation of the IR-safe 
integrals. These packages are surrounded by self-written interface programs, which allow 
for a completely automated computation jllj. As in the first approach, the relevant IR- 
divergent n-point functions are implemented in a Fortran code. They can be extracted 
from Refs. [29 4 5|46] . For the analytical treatment of the abelian five-point functions, 
we refer to Ref. 0Bj- The pentagon diagrams involving a nonabelian coupling can all be 
reduced to integrals with a lower number of external legs by partial fractioning. 

Our analytic result is too lengthy to be presented here. However, the interference of 
the light-quark box amplitude T q with T is sufficiently compact to be listed, which is 
done in Eq. (|A.1|) . 



2.1.1 Parameter and wave- function renormalization 

The self-energy and triangle diagrams of the type indicated in Fig. |21 contain UV sin- 
gularities, which are cancelled upon the renormalizations of the QCD gauge coupling 
g s = y/4:7ia s , the charm-quark mass m and field tp, and the gluon field A^. Specifically, 
the renormalization transformations read 



9° s 



Z g g s 



m 



o 



Z m m, 



Z 3 A 



11) 



where the superscript labels bare quantities and Zi = 1 + 5Zi, with i = g,m, 2,3, are 
renormalization constants. The quantities SZi are of 0(a s ) and contain UV singularities 
and, in general, also finite pieces. The UV singularities are unique, while the finite pieces 
depend on the choice of renormalization scheme. In order to comply with the Lehmann- 
Symanzik-Zimmermann jl^j reduction formula without any additional finite adjustments, 
we define Z2 and Z3 in the on-mass-shell (OS) scheme, the results being 



5Z™ 



5Z? S = — 



-C F ^ 



1 2 , 4vru 2 t 

— + 3 7 £ + 3 In — + 4 



a. 
An 



(/?o - 2C A ) (- 

\ ^uv £ir 



0(e) 



0(e) 



(12) 



where je is Euler's constant and /3q = (11/3)CU — (4/3)Tpn/ is the one-loop coefficient 
of the QCD beta function. We also define m in the OS scheme, by setting 



5Z 



OS 



-3C, 



a s 

47T 



7i; + In — + - + G(e 

euv m 3 



(13) 



However, we adopt the modified minimal-subtraction (MS) scheme to define g s , by putting 

1 



- 7s + ln(47r) 



(14) 
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2.1.2 Operator renormalization 



A crucial feature of effective field theories, such as NRQCD, is that the composite opera- 
tors are generally subject to renormalization. In the case of NRQCD, this is essential in 
order to ensure the complete cancellation of IR and Coulomb singularities at NLO and 
so to overcome the conceptual problems of the CSM mentioned in Section To be con- 
sistent with the rest of our calculation, we also employ dimensional regularization here. 
We adopt the technique described in Refs. |2|31| to directly evaluate the NLO corrections 
to the NRQCD operators. In this way, we avoid having to match partonic cross sections 
evaluated in NRQCD with their counterparts in full QCD. 

In the case under consideration, we have to renormalize the CO ME lO H 3 S^ 



which appears at LO, in \%\ 2 . In d space-time dimensions, this ME has mass dimension 
d — 1. We thus introduce the 't Hooft mass scale of NRQCD, A, to keep its renormalized 
version, which we wish to extract from experimental data, at mass dimension 3. 



The four-quark operator lO H 3 S{ 



is related to the amplitude for the elastic scat- 
tering of a cc pair. The corresponding tree- level diagram is depicted in Fig. EUa). The 
one-loop corrections to this amplitude are obtained by attaching a virtual gluon line in all 
possible ways to the external heavy-quark legs, and they involve self-energy and vertex 
corrections [see Fig. E(b)— (e)]. Using the NRQCD Feynman rules in the quarkonium rest 
frame, expanding the one- loop integrands as Taylor series in 1/m, and performing the 
integration over the loop momentum, we obtain the unrenormalized one-loop result 



O 



H 



35.(8) 



o 



x 



11 



35.(8) 




1 

em 



1 + [Cf - 



Ca 
2 



2 

E 

J=0 



c F (o 



2v 



3p(l) 



+ 



4a s ( Aiifi^ 



37rm 2 



A 2 



exp(-e7£) 



B F (0 



3p(8) 



(15) 



where the subscript labels the tree-level quantity and \x is the 't Hooft mass scale of 
QCD that enters through the <i-dimensional loop integration. The term proportional to 
1/v represents the Coulomb singularity, which arises from the exchange of a longitudinal 
gluon between the outgoing c and c quarks [see Fig. C2c)]. Obviously, NRQCD operators 
of different cc Fock states n start to mix at one loop. Furthermore, the presence of UV 
singularities indicates that they need renormalization. In the following, we choose the MS 
scheme for that. We thus write 



O 



H 



35.(8) 



o 



11 



4a K 



37rm 2 



'47T//V 1 

— — exp{-e-f E ) — 



A 2 



x £ (C F (O h 

J=0 



3p(D 



B F (0 



1-H 



3p(8) 



(16) 



where the subscript r labels the renormalized quantity and we identify A with the NRQCD 
renormalization scale. Inserting Eq. (fTp^l into Eq. (|15p. we obtain 



'O 



H 



''S{ 8) ]) = (0 H [ 3 sJ 8) ]) r (A) 



Ca 
2 



2v 



+ 



4a s I AiTfi 2 



37rm 2 



A 2 



exp(-e7£;)- 



em 
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J=0 



3p(D 



+ B F (0 



(17) 



Combining Eq. (|T7|) with Eq. (JJJ) generates an IR counterterm at 0(a s ) that is indispen- 
sible to render the total NLO result finite. The Coulomb singularity present in Eq. (fTTj) 
is necessary to cancel similar terms in the virtual corrections. 

2.1.3 Renormalization group equation 

The renormalization group equation that determines the A dependence of (O h 3 S[ 8 ^ ^ (A) 
may be derived by differentiating Eq. (fTTj) with respect to A and then taking the physical 
limit e — > 0. It reads 

X 2 d 



dX 2 



o 



11 



35,(8) 



(A) 



4a s (/x) 
37rm 2 



J=0 



3p(l) 



B F (0 



3p(8) 



(18) 



There is no obvious physical reason to distinguish between the scales \i and A, which refer 
to the ccg and cccc vertices in the one-loop diagrams of Figs. E2(b)-(e), respectively. Both 
scales should essentially be of 0(m). In the following, we thus identify fi = A. Integration 
of Eq. (JTj| then yields 



O 



11 



3 S? ] ) (A) 



O 



11 



16 



x f 1 ln a «(^o) 
\/3 a s (A) 



+ 



3m 2 
Pi 



3p(D 



3p(8) 



J=0 



[a«(A) - a s (A )] j 
'35.(8) 



(19) 

where A is a reference scale for which the value of (O h ^Sf^ ^ (A ) is assumed to be 
known and (3\ = (34/3)C^ — ACpTpUf — (2 0/3) C AF F nf is the two-loop coefficient of the 
QCD beta function. In want of a genuine NLO determination of (O h [ 3 S{ 8) ] ) (A ) from 

a fit to experimental data, we choose Ao = m and identify (O h ^ (to) with its 

A-independent LO value, which is known from the literature [2H|. Since we consider the 
NLO corrections to the LO cross section (J7J), which is already of 0(a s ), we have to employ 
in Eq. (|19|) the two-loop formula for a s (X), which reads 

a s (X) 1 (3\Yn.L 



4tt (3 L ffiL 



?r.2 



(20) 



where L 



In (A 2 /A 



2 

QCD 



2.2 Real corrections 

The real corrections to the cross section of process (J2J) arise from the partonic subprocesses 

7 (A;i) + 7 (fc 2 ) -> cc[n](p) + ^(Afc) + 0(fc 4 ), (21) 



10 



where n = 3 Pf\ l sf\ 3 sf , 3 Pf ] , and 

l{ki) + j(fa) -> cc[n] (p) + g (fc 3 ) + g(fc 4 ), (22) 

where n = ^Sq , 3 Si , 3 Pj . The respective diagrams are presented in Figs. Ufa) and EJ 
From the technical point of view, it is convenient to evaluate the gluon polarization sum 
as Spoi £ a(l) £, u(Q) — ~^ ah 9iiv-, at the expense of allowing for the Faddeev-Popov ghosts of 
the gluon, u g and u g , to appear in the final state. The additional partonic subprocesses 
are 

7 (A;i) + 7(/c 2 ) -> cc[n](p) + +%(/c 4 ), (23) 

where n = ^ . The corresponding diagrams are depicted in Fig. ^b). 

Process (j21j) with n = 3 Si is forbidden by Furry's theorem j^H], as may be understood 
by observing that the 3 Si projector effectively closes the charm-quark line and acts like 
a vector coupling and that the two gluons are then in a CS state, so that we are dealing 
with a closed fermion loop containing five vector couplings. This was also verified by 
explicit calculation. On the other hand, process (E2J) with n = 3 S{ 1] , 3 pP is prohibited 
by colour conservation, since the two quark lines are connected by a single gluon, which 
ensures that the cc and qq pairs are both in a CO state. 

The diagrams in Fig. El can be divided into three classes. The diagrams of the first 
(second) class are proportional to e 2 (e 2 ) and contain a timelike virtual gluon that splits 
into a cc (qq) pair. The diagrams of the third class are proportional to e g e c and contain 
a spacelike virtual gluon. 

As for the kinematics of processes (|21| ) - (J23j) . we now have fa + k 2 = p + k 3 + A; 4 , 
k\ — k\ = fc| = k\ = 0, and p 2 = Am 2 . We can define (Vj = 10 Mandelstam variables, 
5 of which are linearly independent. In addition to S\, t\, and u\ defined in Eq. (JHJ), we 
introduce 

S3 = (p+ h) 2 - Am 2 , s 4 = (p+ h) 2 - Am 2 , s 34 = (fc 3 + k 4 ) 2 , 
S13 = (fa - h) 2 , s 2 3 = (fa - fa) 2 , s u = (fa - fa) 2 , s 24 = (fa - fa) 2 . (24) 

If we take Si, t\, U\, s i4 , and s 24 to be the set of independent Mandelstam variables, then 
the others are given by 

S3 = Sl + S U + S24, 

s 4 = — Si — t\ — U\ — s 14 — s 24 — 8m 2 , 
S34 = si + ti + ui + 8m 2 , 
S13 = -si - ti - s u - Am 2 , 

S23 = -s\ -u 1 - s 24 - 4m 2 . (25) 

Generically denoting the T-matrix elements of processes (|21 )l - (j23|l by %, their cross 
sections may be evaluated as 

da r = -^-dPS 3 (fa + fa;p, fa, fa)\%\ 2 . (26) 
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In the case of process (JUJ), a factor of 1/2 has to be included on the right-hand side of 
Eq. ()26|) to account for the fact that there are two identical particles in the final state. 

Integrating Eq. (}26|) over the three-particle phase space while keeping the value of px 
finite, we encounter IR singularities, which can be of the soft and/or collinear type. In 
order to systematically extract these singularities in a Lorentz-invariant way, it is useful 
to slice the phase space by introducing infinitesimal dimensionless cut-off parameters 5i 
and 5f, which are connected with the initial and final states, respectively |48j . In the case 
of processes (J2Tj) and we are thus led to distinguish between the following regions of 
phase space. In the soft region, where either S3 + S34 < 5fS or s 4 + S34 < 5fS, one of the 
outgoing gluons is soft and, in addition, they may be collinear. In the final-state collinear 
region, where s 34 < 5fS < s 3 +s 34l , S4+S34, the outgoing gluons are both hard and collinear. 
In the hard region, which comprises the residual phase space, none of the outgoing gluons 
is soft and they are not collinear. In the case of process ()22|). we differentiate between 
the following regions of phase space. In the initial-state collinear region, where one of the 
inequalities S13 < <5jS, s u < 5iS, s 23 < Si,s, and s 24: < <5,s is satisfied, one of the outgoing 
q and q quarks is collinear to one of the incoming photons. In the final-state collinear 
region, where s 34 < 5js < s 3 + s 34 , s 4 + s 34 , the outgoing q and q quarks are both hard 
and collinear. In the hard region, which includes the left-over phase space, none of the 
outgoing q and q quarks is soft and they are not collinear. While the various contributions 
depend on the cut-off parameters 5i and/or 5f, their sum must be independent of them. 
The verification of the cut-off independence provides a nontrivial and powerful check for 
the correctness of our calculation. 



2.2.1 Hard region 

The integration of Eq. over the hard region can be facilitated by decomposing the 
three-particle phase space as 

dPS 3 (k 1 + k 2 ;p,k 3 ,k i ) = — dPS^&i + k 2 ; p, k u )dPS 2 (k M ; k 3 , fc 4 ), (27) 

Z7T 



where k 34 = k 3 + k 4 and 

A-d d d P x , 2 a^2\q/ 0\ <^34 



dPSt{h + k 2 - p, k u ) = ^—J^Sip 2 - Am 2 )d{p )^-p [ {27i) d 5^ d \k 1 + k 2 -p- k 3A ) 

1 / 4ttu 2 s Y , , 

= ^7 r , o dhdut. 28 

87rsr(l - e) \t lUl - 4m 2 s J 

On the other hand, writing 

p"=(p°,5r,p), 

k 3 = (k®, . . . , &3 sin 6*3 sin 3 , k 3 sin ^3 cos 03, k 3 cos 9 3 ^j , 

k% = (h®, . . . , &4 sin 6*4 sin 4 , sin # 4 cos 4 , k® cos 64J , (29) 
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we have 

T(l-e) /47r/! 2 Y s + 4m 2 - 2y/sp° 



dPS 2 (k 34 ; k 3 , fc 4 ) 



167r 2 r(l-2e) V(^) 2 / (Vs -P° + 2pcos6 3 ) 2 

x sin d - 3 6 3 d6 3 sin^" 4 3 d0 3 , (30) 



where the boundaries of integration are fixed by the kinematical constraints specified 
above. Since the contribution from the hard region is UV and IR finite, we may perform 
the phase-space integration in d = 4 space-time dimensions. 

2.2.2 Soft region 

As an example, let us consider the kinematical situation where k 3 in processes (|2ip or 
(I23|) becomes sufficiently soft, so that it does not influence the kinematics of the hard 
process involving the residual four-momenta. In practice, this is arranged by imposing 
the condition s 3 + s 34 < 5fS with a sufficiently small value of 5f. The three-particle phase 
space can then be decomposed as 

dPS 3 (h + k 2 ; p, k 3 , h) = dPS 2 {k 1 + k 2 ; p, fc 4 )cffS 3 , (31) 

where 

Jd— 1 I, 

dPS 3 s = ^- d d k3 



(2n) d - 1 2k^ 
r(l - e) / rtfi 2 



87r 3 r(l - 2e) 



kldkl sin d - 3 3 d6 3 sin d - 4 <p 3 d<p 3 . (32) 



At the same time, IXI 2 factorizes as 



|T r | 2 = |T | 2 F 3 , (33) 

where F 3 is an appropriate Eikonal factor, the specific form of which depends on the 
cc Fock state n. Specifically, soft singularities occur for processes (}2"Tj) and (J2l?j) with 
n = 3 Pj^ , 3 S[ 8 ^ , 3 Pj 8 ' ) • The appropriate expressions for F 3 are listed in Table El Notice 
that, in Eq. (j^H)), T always refers to process (j2J) with n = 3 Si . 

The integration of the Eikonal factors over the soft parts of the phase space can be 
performed analytically, and the results for the soft integrals, 



Ij= dP$lF 3 (34) 

Jk°<8 f v r s/2 



and similarly for fc 4 , may be found in Table El where 



P= S ~^r (35) 
s + 4m 2 v ' 
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cess 



77 
77 

77 

77 
77 



cc 



cc 



cc 



cc 



apPT 



(8) 



3 S[ 



Sp(8) 



3 S[ 



(8) 



qq 



cc[n\qq 



Limit 



h 
h 

k 4 

h 
h 

k 4 

h 

k 2 
k 2 



-> 
-> 
-> 
-> 
fc 4 
-> 
-> 
fc 4 
A; 3 

fc 4 



Eikonal factor 



^3 = 

F 4 = 
F 3 = 
F 4 = 
F34 = 
F 3 = 
F 4 = 
F34 = 
F23 = 
F24 = 



53 V s 34 

54 \_S34 

«34 
16B F g 2 s 
3si 



4m 
S3 
4m 2 
S4 
1-z 



3 2 ( 1 -f) 

Wz& (l _ |) 



+ (1 - - 2ez(l - z)\ 



«34 L 

_2e|e2 2 2 +(1 _ 2) 2_ e 



S?4 2 



Table 2: Eikonal factors appropriate for the soft and collinear limits of the squared T- 
matrix elements of the various partonic subprocesses. Here, n = ^q 8 ^, 3 s[ 8 \ 3 Pj . 



2.2.3 Final-state collinear region 

In the final-state collinear region, where S34 < SfS < s 3 + s 34 , S4 + S34, it is useful to 
introduce the parameterization 

K = ((1 - z)p + 2{1 ^ z y -kr, -(1 - z)p ] j , (36) 

where z, with < z < 1, determines the fractions of longitudinal momentum that the final- 
state partons receive from the splitting gluon, k\ = z{\ — -2)534, and terms of O (k^/p 4 ) 
are neglected. We can then decompose the three-particle phase space as 

dPS 3 (h + k 2 ;p, k 3 , fc 4 ) = dPS 2 (k 1 + k 2 - p, k u )dPSf, (37) 

where 

^34 _ ,,4-d d d 1 /C3 /C34 



^/p034 _ A 

bc (2^2^ k\ 

1 ' { dzds 34: . (38) 



167r 2 r(l-e) V^(l- ^)s 3 4 

Using the relation 



r(l-e) r(l-2e) 
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L ' (1 +0{e% (39) 



Process 



Soft and collinear integrals 



77 — > cc 



77 — > cc 



77 — ► cc 



77 — ► cc 



3p(D 



3^(8) 



3p(8) 



3g(8) 



99 



99 



99 



99 



J 3 

s 

J 4 

a 



3m 2 



2 111 i-fl ^ a 1 1 



If = -2C A g 2 C e (0 [i (21n^ 



/ 3 

s 

i; 1 



8Bf£ c 



3m 2 



+ ln 2 ^ 

<5 /s 



f I + I i n 2±£ _ I 

V e ^ 8 111 1-8 3 



^ = -\T t n f9 tC^{\ + \) 



67 
18 



Table 3: Soft and collinear integrals. 



where the terms of 0(e 2 ) are inconsequential in the physical limit d — > 4 because, in the 
collinear limit, \%\ 2 only develops simple poles in e, we recover the familiar combination 
of Euler's Gamma functions. At the same time, 171 1 2 factorizes as 



\%\ 



\%\ 2 F, 



3 b 



(40) 

where F34 is an appropriate eikonal factor. In the case under consideration, final-state 
collinear singularities occur for processes (J2l|) and (}22|) with n = 3 s[ 8 \ where they are 



related to g — > gg and g — > qq splitting, respectively. The corresponding expressions for 
F34 may be found in Table 121 They are related to the timelike g — > g and g — > q splitting 
functions, 

l-z 



as 



-P34 = 
-P34 = 



2C A 



1 



+ 



+ z(l-z) 



z 2 + {l-zf 



(41) 
(42) 



S34 

M [Pj ^( z ) _ 2eT F z(l - *)] 
S34 



(43) 



respectively. 

The integration of F34 over the collinear parts of the phase space can be performed 
analytically, and the resulting collinear integrals, 



r 34 



S34<Sfi 



rfPSfF 34 , 



(44) 



may be found in Table El In the case of process (|21|). the range of integration in Eq. 
has to be further constrained by the condition Sfs/si < z < 1 —Sfs/si to ensure that the 
final-state gluons are both hard. 
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2.2.4 Initial-state collinear region 

As an example, let us consider the kinematical situation where the incoming photon with 
four-momentum k 2 in process f)22j) splits into a hard, collinear, quasireal antiquark q, with 
four- momentum k' 2 , and the outgoing quark q, with four-momentum ks. The collinearity 
is enforced by the constraint s 2 3 < 5iS with a sufficiently small value of <5j. It is then 
helpful to introduce the parameterization 



K = (k 2 , T , k 2 



k'2 = {^zk 2 + j k T , zk 2 ^ , 

K = ((1 - z)k 2 + 2{1 k } z)h , ~k T , (1 - z)k 2 j , (45) 

where k\ = —(1 — z)s23 an d terms of O (k T /k 2 ) are neglected. We can then decompose 
the three-particle phase space as 

dPS 3 (h + k 2 ; p, k 3 , h) = dPS^h + zk 2 ; p, k 4 )dPS 23 , (46) 

where 

Ad- It. 



(27r) Q! - 1 2fc§ 
1 / 4vr/i 2 



16vr 2 r(l-e) \-(l-z)s 23 



dzds 2 %. (47) 



At the same time, \7~ r \ 2 factorizes as 



|T r |2 = |T |2 (7^) +q(k' 2 ) -> ds[n](p) +g(A; 4 ))F 23 , (48) 

where the cc Fock state n is the same as appears in % and F 23 is an appropriate Eikonal 
factor. In the case under consideration, initial-state collinear singularities occur for pro- 
cess (J2*2*|) with n = ^iSq , , ^Pj , where they are always related to 7 — > qq splitting. 
The corresponding expressions for the Eikonal factors Fij, where i = 1,2 and j = 3,4, 
may be found in Table El They are related to the spacelike 7 — > q splitting function, 

P^ q {z) = ^—P g ^ q (z), (49) 

J-F 

where P g -+ q (z) is defined in Eq. (|4*2*j) . as 



2e 2 

F . = - — 

ZSij 



P 7 _> 9 (Z) ^ 2 

iV, 



(50) 
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The integration of Eq. (jl8|) over the collinear region of the phase space, with s 2 3 < SiS, 
can be performed analytically. In terms of differential cross sections, the result reads 

47T/! 2 X 6 



da 



23 




- z)8 i s j 

x da(-i(k l ) + q(zk 2 ) -* cc[n}(p) + q(k 4 )), 



[R^(z) - ee iV c ] 



(51) 



where s is the two-photon invariant mass square defined above Eq. (jHJ). According to the 
mass factorization theorem j3H], the form of this collinear singularity connected with an 
incoming-photon leg is universal and can be absorbed into the bare PDF of the antiquark 
q inside the resolved photon, fq/j(z). To this end, one defines the renormalized photon 
PDF as 



a 



1 rfi 



r l-2e 



P^ q {z) + ep^z) 



(52) 



where M is the factorization scale and the finite function CL 



is introduced to fix 



the factorization scheme. In the modified minimal-subtraction (MS) scheme, we have 

in Ref. [15] . we have 



Cj(z) = 0, while in the so-called DIS 7 scheme, which is used, e.g 



Cy(z) = N c 



z 2 + 



l-z) 



. l-z 
In h 



\z(l-z) 



1 



(53) 



Incorporating the second term on the right-hand side of Eq. ()52|) in Eq. (JSTj) . we are left 
with the finite contribution, 



da 



23 



dz 



a: 
2^ 



Pj^q(z) In • 



z)SiS 



+ e 2 q [N c - C 7 (z)} 



x da(-i(k l ) + q(zk 2 ) 



M 2 

cc[n}(p) +q{k 4 )) 



(54) 



which depends on the choice of M. In turn, the resolved-photon contribution is evaluated 
with the renormalized photon PDF fq/~,(z, M 2 ), which also depends on M. In the sum 
of these two contributions, this M dependence cancels up to terms beyond NLO. In this 
sense, the notion of direct-photon contribution ceases to be separately well defined at 
NLO. 



2.3 Assembly of the NLO cross section 

The NLO result for the cross section of process is obtained by adding the virtual 
and real corrections to the LO result. Collecting the various contributions discussed 
above, arising from the parameter and wave-function renormalization (ct), the operator 
redefinition (op), the initial-state (is) and final-state (fs) collinear configurations, the soft- 
gluon radiation (so), and the hard-parton emission (ha), we can schematically write the 
resulting differential cross section of 77 —>■ H + X, including the MEs, as 

da(n, A, M) = da (n, A)[l + 5 vi (fi; e w , em, v) + 5 ct (/x; euv, em) + <$op(^, A; e m , v) 
+ 6fs0"; em, 5 f )} + da is (n, A, M; 8^ + da so (fi, A; e IR , S f ) 
+ da ha (p:,X;Si,Sf), (55) 
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where the dependences on the unphysical mass scales /i, A, and M and the regulators era, 
euv; v i an d Sf are indicated in parentheses for each term. Everywhere in Eq. (J55j) . 
a s {jj) is evaluated in the MS scheme using the two-loop formula (J20J), m is defined in 
the OS scheme, and the MEs are understood as their MS values (0 H [n]) r (X). This is 
necessary to ensure the exact cancellation of the singularities. 

The right-hand side of Eq. ()55|) is manifestly finite. The UV divergences cancel be- 
tween <5 vi and <5 ct ; the IR singularities among <5 vi , S ct , S op , Sf s , and da so ; and the Coulomb 
singularities between 5 V ; and S op . Note that S op is UV finite upon operator renormalization 
and that d<j m is IR finite upon photon PDF renormalization; therefore, 5 op and da- ls in 
Eq. (jo3|) do not depend any more on euv an d £ir, respectively. The full details about 
the cancellation of the IR singularities in the various cc Fock states n are presented in 
Table EJ where the following short-hand notation is used: 



Subprocess 



3^W 



Source 



IR-singular term 



77 

77 

77 

77 

77 
77 
7<? 



ex- 



ec 



cc 



ccfpf] + gg 



cc 



+ 9 



+ 99 



3 S[ 8) ] + 



99 



cc[n] + qq 
cc[n] + q 



qq 



virtual 

operator 
operator 
soft 

final-state 
soft 

operator 
soft 

final-state 
initial-state 
mass fact. 



+ 



4tt^ 
m? 

3e 



exp(-e7£;) 



2e 111 2m 2 + 12 



|T n | 2 



3^(8) 



3p(D 



f^(¥) e - P (-67.)H^ 
^(^) e exp(- e7 ,)l(i-lnfe + l)|T | 

-^(#) £ -p(-e 7 .)i|T p 



3p(8) 



( (i-z)8 iS ) f .,,(•:)'/'-„ /' 



1 Af 2 



±p 



7" 



7^9 
+5 (^) T^q^ccq 



■yq^ccq [ 

n] 



Table 4: Compilation of the IR-sinj 
various partonic subprocesses. Here, 



;ular terms arising from the various sources in the 

„ _ lc(8) 3o(8) 3p(8) 



The right-hand side of Eq. (J55j) is independent of the cut-off parameters Si and Sf. 
Specifically, the Si dependence cancels between da- m and da^ and the Sf dependence among 
Sf s , da so , and da^. While Sf s , da lS) and da SQ are known in analytic form, the phase-space 
integrals occurring in the evaluation of da^ are rather cumbersome and are thus solved 
numerically. Consequently, the cancellation of Si and Sf has to be established numerically, 
too. The goodness of this cancellation is assessed in Section El 
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The right-hand side of Eq. is independent of the renormalization scales /i and 
A up to terms that are formally beyond NLO; this cancellation is not exact because 
the running of a s {fi) and (0 H [n]) r (X) is determined from the respective renormalization 
group equations, which resum logarithmic corrections to all orders. However, the right- 
hand side of Eq. depends on the factorization scale M already in NLO; this M 
dependence is cancelled up to terms beyond NLO once the cross section of single-resolved 
photoproduction, which depends on M through the photon PDFs f a / 7 (x a , M), is added. 
If M is to be varied, this must thus be done simultaneously in the cross sections of direct 
and single-resolved photoproduction. 



3 Numerical results 

We are now in a position to present our numerical analysis of the inclusive production 
of prompt J/ip mesons in two-photon collisions with direct photon interactions at NLO 
in the NRQCD factorization framework. We consider TESLA in its e + e~ mode with 
y/s = 500 GeV, where the photons are produced via bremsstrahlung and beamstrahlung. 

The discussion of the numerical analysis proceeds in three steps. We first specify our 
input parameters. We then verify that our numerical evaluation is independent of the 
technical cut-off parameters 5{ and 5f and assess its dependence on the renormalization 
and factorization scales fi, A, and M. Finally, we explore the phenomenological conse- 
quences of our analysis by studying the size and impact of the NLO corrections on the 
Pt and y distributions of the cross section. 



3.1 Input parameters 

We use m = 1.5 GeV and a = 1/137.036. For direct photoproduction at NLO (LO), we 
employ the two-loop (one-loop) formula for a^ f \fi) with rif = 3 active quark flavours 
and Aqq D = 299 MeV (204 MeV) gSj. We use the photon PDFs from Gliick, Reya, and 
Schienbein (GRS) jlHI, which are the only available ones that are implemented in the 
fixed-flavour-number scheme, with rif = 3. When we combine the contribution due to 
single-resolved photoproduction, which is so far only known at LO, with the NLO one due 
to direct photoproduction, then we nevertheless evaluate it using the NLO formula for 
ai lf \fi) and the NLO set of the GRS photon PDFs so as to maximize the compensation 
of the M dependence. In order to render this choice fully consistent, we will have to in- 
clude the NLO correction to the single-resolved contribution and also the double-resolved 
contribution at NLO, once they become available. On the other hand, we consistently 
evaluate the single-resolved contribution using the LO formula for ai ll \fi) and the LO 
set of the GRS photon PDFs when we consider it separately or in combination with 
the direct contribution at LO. Our default choice of renormalization and factorization 
scales is fi = M = and A = m. In want of NLO sets of J/ip, Xcj, and ip' MEs, we 
adopt the LO sets determined in Ref. using the LO set of proton PDFs from Mar- 



tin, Roberts, Stirling, and Thorne (MRST98LO) jSO]. Specifically, (C^(nS) 3 S{ 



and 
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3d(1) 



and Xc2 
and / O Xc0 



were extracted from the measured partial decay widths of ip(nS) — ► / + Z 



77 



respectively, while (0^ nS ^ 



15.(8) 



Qi>(nS) 3gfj- 



t(8) 



3p(8) 



3 5'| 8 ' 1 ) were fitted to the transverse-momentum distributions of ip(nS) and x 



cJ 



inclusive hadropro duct ion [5| and the cross-section ratio cr Xc2 /o" Xcl [SI] measured at the 



Tevatron. The fit results for (p^ nS) 
so that the linear combination 



15.(8) 



and ( 0^ nS ^ 



3p(8) 



are strongly correlated, 



M i>(nS) 



_| / Q^(jiS) 

m? ^ 



3 P, 



(8) 



(57) 



with a suitable value of r, is quoted. Unfortunately, Eq. ()55|) is sensitive to linear com- 



bination of 
in Eq. (IB7I) 



£)V(™S) 



1 SP ) and 



In want of more specific in: 



3 P f 



(8) 



Qil>(nS) 



15.(8) 



(r/m 2 ) 



3p(8) 



that is different from the one appearing 
brmation, we thus make the democratic choice 

We now discuss the photon flux functions that enter our predictions for photoproduc- 
tion in the e + e~ mode of TESLA. The energy spectrum of the bremsstrahlung photons 
is well described in the WWA by Eq. (27) of Ref. [52]. We assume that the scattered 
electrons and positrons will be antitagged, as was usually the case at LEP2, and take the 
maximum scattering angle to be # max = 25 mrad jSS]. The energy spectrum of the beam- 
strahlung photons is approximately described by Eq. (2.14) of Ref. [54J. It is controlled 
by the effective beamstrahlung parameter T, which is given by Eq. (2.10) of that refer- 
ence. Inserting the relevant TESLA parameters for the y/S = 500 GeV baseline design 
specified in Table 1.3.1 of Ref. [53] in that formula, we obtain T = 0.053. We coherently 
superimpose the WWA and beamstrahlung spectra. 

With the NLO corrections to single- and double-resolved photoproduction yet to be 
evaluated, we are not in a position to present a complete phenomenological prediction 
that could be confronted with experimental data as it stands. Therefore, we refrain from 
presenting a full-fledged quantitative estimate of the theoretical uncertainties. However, 
in the next section, we do investigate the dependences on the renormalization and factor- 
ization scales /1, A, and M. 



3.2 Academic study 

For the purpose of the following technical study, it is sufficient to consider a typical 
kinematic situation. We thus choose as our reference quantity the differential cross section 
d 2 cr/dp T dy at px = 5 GeV and y = 0. 

In Figs. EJa) and (b), the NLO result of direct photoproduction (solid lines), its hard, 
noncollinear component (dot-dashed lines), which corresponds to da^ in Eq. (|55jh and 
the remainder (dotted lines) are plotted as functions of 5i for 5f = 5 x 10~ 3 and as 
functions of 5f for 5i = 5 x 10~ 4 , respectively. For comparison, also the LO result (dashed 
lines), which is, of course, independent of 5i and 6j, is shown. Notice that Si and Sf 
are varied over several orders of magnitude. They have to be chosen judiciously: if they 
are chosen too large, then the collinear approximation underlying Sf s and d<Ji S and the 
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soft approximation underlying da so break down; if they are chosen too small, then the 
numerical evaluation of da^ becomes inaccurate. Our default values for the rest of the 
numerical analysis are 5i = 5 x 10~ 4 and 5f = 5 x 10~ 3 . Since the remainder depends 
logarithmically on 5i and 5f, it is represented by straight lines in Figs. EUa) and (b), the 
abscissae of which are logarithmic. The same is true, to good approximation, for the hard, 
noncollinear component, and the combined NLO result is practically independent of 5i 
and Sf, as it should. This suggests that the numerical integration is rather precise and 
stable. We conclude that this error source can be safely neglected in the determination 
of the theoretical uncertainty of the NLO result. 

In Fig. the NLO (solid line) and LO (dashed line) contributions due to the cc 
Fock state n = ^j 8 '* are shown as functions of fi, while A and M are kept fixed at their 
reference values. We focus on this particular channel because it is the only one already 
open at LO, so that a compensation of the \x dependence can occur. Notice that fi is 
varied by more than one order of magnitude, from m-r/4 to Arrir- Passing from LO to 
NLO, the \i dependence is appreciably reduced, reflecting the partial compensation of 
the \i dependence of da by the one of 5 ct in Eq. ()55|). In fact, the related theoretical 
uncertainty amounts to ^3q% at LO and to ±18% at NLO. 

In Fig. |HJ the NLO result is shown as a function of A (solid line), while \i and M 
are kept fixed at their reference values. Notice that A is varied by almost one order 
of magnitude, from m/2 to 4m. The related theoretical uncertainty amounts to ±50%. 
The A-dependent terms in Eq. (|53|) are dao and 5 op . In order to exhibit the partial 
compensation in A dependence between these two terms, we also include in Fig. |H] the 
results that are obtained by only varying A in da^ (dashed line) or 5 op (dotted line) at 
a time. The dotted line is straight, reflecting the fact that 5 op depends on A through a 
single logarithm. 

In Fig. IH1 the NLO result of direct photoproduction (dotted line), the LO result of 
single-resolved photoproduction evaluated with the NLO versions of aia ((jl) and the 
photon PDFs (dashed line), and their sum (solid line) are shown as functions of M, 
while \x and A are kept fixed at their reference values. Notice that M is varied by more 
than one order of magnitude, from m^/4 to 4m^. The M dependences of the individual 
contributions almost cancel each other, leaving a theoretical uncertainty of ±11% on their 
sum. 

3.3 Phenomenological study 

In Fig. HI3 we study d 2 a/ dpr dy (a) for y = as a function of pt and (b) for pr = 5 GeV as 
a function of y. In each case, the LO (dashed line) and NLO (solid line) results of direct 
photoproduction as well as the LO result of single-resolved photoproduction evaluated 
with the LO versions of a^ f \fj.) and the photon PDFs (dotted line) are shown. Notice 
that our analysis is only valid for finite values of pt', in the limit px —> 0, additional 
IR singularities occur, which require a more sophisticated scheme of phase space slicing. 
Therefore, we do not consider p T values below 2 GeV in Fig. [TOT a). For a more detailed 
discussion of this point, we refer to Ref. [2*T] . 
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From Fig. 110( a). we observe that, with increasing value of px, the NLO result of direct 
photoproduction falls of considerably more slowly than the LO one. This is also evident 
from Fig. ITlT a). where the NLO to LO ratio is shown for y = as a function of px 
(solid line). This feature may be understood by observing that so-called fragmentation- 
prone partonic subprocesses (SEj start to contribute to direct photoproduction at NLO, 
while they are absent at LO. Such subprocesses contain a gluon with small virtuality, 
q 2 = Am 2 , that splits into a cc pair in the Fock state n = 3 S[ 8 ^ and thus generally generate 
dominant contributions at px 3> 2m due to the presence of a large gluon propagator. In 
the case under consideration, the relevant T-matrix elements are those of 77 — > cc 3 s[ 8 ^ g 

and 77 — > cc 3 S[ 8 ^ qq that are proportional to e 2 ; the corresponding diagrams are the 
last two ones in Fig. |21 and the last six ones in Fig. El respectively. In single-resolved 
photoproduction, a fragmentation-prone partonic subprocess already contributes at LO; 
the relevant T-matrix elements is the one of ^q — > cc 3 S[ 8S> q that is proportional to e q . 
This explains why the solid and dotted curves in Fig. HOT a) run parallel in the upper p T 
range. At low values of px, the fragmentation-prone partonic subprocesses do not matter, 
and the relative suppression of direct photoproduction is due to the fact that, at LO, this 
is a pure CO process. 

At pt = 5 GeV, single-resolved photoproduction is still overwhelming, as is evident 
from Fig. HOT b). In Fig. fTOT b). the two pronounced maxima in the dotted line may be 
traced to the T-matrix element of 75 — > cc[n]g that is proportional to e c , which contains 
a virtual gluon in the t channel that can become almost collinear with the incoming q 
quark, and to the one that is proportional to e q , which contains a virtual q quark in the 
u channel that can become almost collinear with the incoming photon. 

Furthermore, we observe from Fig. ITTT b) that the NLO correction to direct photopro- 
duction dramatically increases towards the forward and backward directions. This feature 
is also nicely exhibited in Fig. ITTT b). where the NLO to LO ratio is shown for px = 5 GeV 
as a function of y (solid line). This is due to the finite remainders of the initial-state 
collinear singularities that were absorbed into the photon PDFs. 

In Fig. HH two kinds of QCD correction (K) factor are studied (a) for y = as a 
function of px and (b) for px = 5 GeV as a function of y. The one that is defined as 
the ratio of the NLO and LO results of direct photoproduction (solid lines) was already 
discussed above. An alternative definition of K factor is obtained by adding the LO re- 
sult of single-resolved photoproduction to both numerator and denominator, where it is 
evaluated with the NLO and LO versions of aj 1 (fx) and the photon PDFs, respectively 
(dashed lines). The rationale for this definition is the circumstance that the NLO result 
of direct photoproduction and the LO result of single-resolved photoproduction are in- 
terconnected by mass factorization as explained in Section 03.4. For the reasons given 
in Section Ell, this definition of K factor is slightly inconsistent and the resulting values 
should, therefore, be taken with a grain of salt. Since the LO result of single-resolved 
photoproduction dominates both numerator and denominator of this K factor, the latter 
takes moderate values. The fact that these values are even below unity in the considered 
regions of phase space is due to the different choices of a£ nj (fj) and the photon PDFs for 
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the single-resolved contribution in the numerator and denominator of this K factor. 

4 Conclusions 

In the case of charmonium, the experimental verification of the NRQCD factorization 
hypothesis in two-photon collisions at LEP2, photoproduction at HERA, and hadropro- 
duction at the Tevatron has now practically come to a halt, one reason being that the 
available theoretical predictions for the production of prompt J / ip mesons with finite val- 
ues of Pt in high-energetic photon-photon, photon-hadron, and hadron-hadron collisions 
are only of LO and thus suffer from considerable uncertainties, mostly from the depen- 
dences on the renormalization and factorization scales and from the lack of information 
on the nonperturbative MEs. 

Having recently considered the cross sections of processes (|2]) [21] and (}2Tf - ([23|) [T9] 
in connection with the associated production of prompt J/ip mesons with one or two 
identified energetic hadron jets, respectively, in this paper, we took the next step by 
studying at NLO the inclusive production of prompt J/ip mesons with finite values of 
Pt- This is the first time that an inclusive 2 — > 2 process was treated at NLO in the 
NRQCD factorization framework. The technical difficulties that needed to be tackled 
include the treatment of UV, collinear, soft, and Coulomb singularities, the NRQCD 
operator renormalization, the mass factorization, and the analytic evaluation of five-point 
one- loop integrals with UV, soft, and Coulomb singularities. 

As for the real corrections, we employed the phase-space slicing method to demarcate 
the regions of phase space containing soft and collinear singularities from the hard regions, 
where the phase-space integrations were carried out numerically. We verified that the 
combined result is, to very good approximation, independent of the choices of the cut- 
off parameters 5i and 5f, over an extended range of values. We worked in dimensional 
regularization in connection with the MS renormalization and factorization schemes, so 
that our NLO result depends on the QCD and NRQCD renormalization scales fi and A, 
respectively, and on the factorization scale M connected with the collinear splitting of the 
incoming photons into massless qq pairs. While the /i and A dependences are formally 
cancelled up to terms beyond NLO within direct photoproduction, the M dependence 
is only compensated by the LO cross section of single-resolved photoproduction. By 
the same token, the strong M dependence of the latter is considerably reduced by the 
inclusion of our new result. This is a crucial phenomenological merit of our work, as may 
be appreciated by observing that, at LO, the overwhelming bulk of the cross section of 
prompt J I ip production in two-photon collisions is due to single- resolved photoproduction 

nun]. 

The K factor of direct photoproduction turned out to be very substantial at large val- 
ues of pt because fragmentation-prone channels start to open up at NLO [see Fig. II 1( a)] . 
In fact, at pt 3> 2m, the pt distribution of direct photoproduction at NLO is rather sim- 
ilar to the LO result of single-resolved photoproduction both in shape and normalization 
[see Fig. HUta)]. 
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In order to complete the NLO treatment of prompt J/ip production in two-photon 
collisions, we still need to evaluate the NLO corrections to single- and double-resolved 
photoproduction. Then, also prompt J/ip production in photoproduction at HERA and 
hadroproduction at the Tevatron can be described at NLO. This will provide a solid basis 
for an ultimate test of the NRQCD factorization framework. 
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A Light-quark loop contribution 



The interference of the light-quark box amplitude T q with the tree-level one %, which 
enters Eq. (JXDJ), reads 



2Re(T *T g ; 



16ey c eY s (0[ 3 S[ 8) [ 
3m 



' stu[s 2 + t 2 + 3tu + u 2 + 3s(t + u)] 
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128m 4 
(s + tfis + ufit + uf 
+ t 2 u 2 (t 2 + u 2 ) + s 2 (t 4 - At s u - Ut 2 u 2 - Atu 3 + u 4 )] J B (4m 2 , 0, 0) 

+ «)+.(«■ + «')] 
(s + t)(s+u)(t + uf <"->'> 
+ «[At-4»)-W + «.1 _ 



+ 



lu[-Ast 2 + t 2 u + s 2 (-4t + u)} 
(s + t) 3 {s + u)(t + u) 



5 (u,0,0) 



4s(2s 2 + t 2 + u 2 ) „ , 
+ , ~ ~ - ~ r C (0, 0, s, 0, 0, 0) 



+ 



(s + t)(s + u)(t + u) 
At(s 2 + 2t 2 + u 2 ) 

(s + t)(s + u)(t + u) 
Au(s 2 + t 2 + 2u 2 ) 

(s + t)(s + u)(t + u) 

A(s 2 + 2t 2 + u 2 ^ 



C (0,t,0,0,0,0) 
Q,(0, 0,0, 0,0) 



(s + t)(t + u) 
A{s 2 + t 2 + 2u 2 ) 
(s + u)(t + u) 



C (4m 2 ,0,t, 0,0,0) 
C (Am 2 ,0,u, 0,0,0) 



24 



4(2s 2 + t 2 + u 2 
(s + t)(s + u) 
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(A.l) 



where B , C , and D denote the one-loop scalar two-, three-, and four-point functions 
[57] in the notation of Ref. j^Hj (see also Ref. [515]). Since the C and D functions are IR 
divergent, we list them in analytical form. We have 

B {s, 0,0) = C t (I-lnZf + 2 



C (s, 0,0, 0,0,0) = C- 
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si 

+ 2 Li 2 I 1 
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e 4m 4 
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4m 4 



Am 2 
3C(2) 



(A.2) 



where £(2) = 7r 2 /6, s = s + ze, i = t + ie, and C e is defined in Eq. The result for 

C (s, 0, 0, 0, 0, 0) can be found, e.g., in Ref. jlHj, and those for C (4m 2 , 0, s, 0, 0, 0) and 
_D (4m 2 , 0, 0, 0, s, t, 0, 0, 0, 0) can be extracted from Ref. 
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Figure 1: Tree- level Feynman diagrams pertinent to the partonic subprocess (j2J). 
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(a) 

Figure 6: Differential cross section d 2 a/dpT dy in pb/GeV of e + e~ — > e + e~J/ip + X in 
direct photoproduction at TESLA with y/s = 500 GeV for prompt J/ip mesons with 
Pt = 5 GeV and y — 0. The NLO result (solid lines), its hard, noncollinear component 
(dot-dashed lines), which corresponds to da^ in Eq. (J33J), and the remainder (dotted 
lines) are shown (a) as functions of Si for Sf = 5 x 10~ 3 and (b) as functions of Sf for 
Si = 10~ 4 . For comparison, also the LO result (dashed lines) is shown. 
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Fig. El (continued). 
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Figure 7: Differential cross section d 2 a/dp T dy in pb/GeV of e + e~ — > e + e^ J/ip + X in 
direct photoproduction at TESLA with y^s = 500 GeV for prompt J ftp mesons with 
p T = 5 GeV and y — 0. The NLO (solid line) and LO (dashed line) contributions due to 
the cc Fock state n = 3 S[ 8 ^ are shown as functions of fi. 
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Figure 8: Differential cross section cPa/dpT dy in pb/GeV of e + e~ — > e + e~ J/ip + X in 
direct photoproduction at TESLA with y/s = 500 GeV for prompt J/ip mesons with 
Pt = 5 GeV and y = 0. The NLO result is shown as a function of A, (i) keeping A in 
the partonic cross sections fixed (dashed line), (ii) keeping A in ((D H 3 S[ 8 ^ ^ (A) fixed 
(dotted line), and (iii) varying all occurrences of A simultaneously (solid line). 
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Figure 9: Differential cross section cPa/dpr dy in pb/GeV of e + e~ — > e + e~ J/ip + X at 
TESLA with ^/s = 500 GeV for prompt J/ip mesons with p T = 5 GeV and y = 0. 
The NLO result of direct photoproduction (dotted line), the LO result of single-resolved 
photoproduction (dashed line), and their sum (solid line) are shown as functions of M. 
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Figure 10: Differential cross section d 2 a / dpr dy in pb/GeV of e + e~ — > e + e~ J/ip + X at 
TESLA with = 500 GeV for prompt J ftp mesons (a) with y = as a function of 
Pt and (b) with px = 5 GeV as a function y. The LO (dashed line) and NLO (solid 
line) results of direct photoproduction are compared with the LO result of single-resolved 
photoproduction (dotted line). 
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Fig. [TUl (continued). 
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Figure 11: QCD correction factor K of the differential cross section cPa/dpr dy of e + e~ — > 
e + e-J/ip + X at TESLA with v/s = 500 GeV for prompt J/V> mesons (a) with y = 
as a function of p T and (b) with p T = 5 GeV as a function K is defined as (i) the 
ratio of the NLO and LO results of direct photoproduction (solid lines) and (ii) the same 
quantity after adding the LO result of single-resolved photoproduction to both numerator 
and denominator, where it is evaluated with the NLO and LO versions of ai n/ ^(/x) and 
the photon PDFs, respectively (dashed lines). 
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